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In this paper, we present charged dilatonic black holes in gravity’s rainbow. We study geometric 
and thermodynamic properties of black hole solutions. We also investigate the effects of rainbow 
functions on different thermodynamic quantities for these charged black holes in dilatonic gravity’s 
rainbow. Then, we demonstrate that hrst law of thermodynamics is valid for these solutions. After 
that, we investigate thermal stability of the solutions using canonical ensemble and analyze the 
effects of different rainbow functions on thermal stability. In addition, we present some arguments 
regarding the bound and phase transition points in context of geometrical thermodynamics. We 
also study the phase transition in extended phase space in which cosmological constant is treated as 
the thermodynamic pressure. Finally, we use another approach to calculate and demonstrate that 
obtained critical points in extended phase space are representing a second order phase transition for 
these black holes. 


I. INTRODUCTION 


Motivated by works on Lifshitz scaling in condensed matter physics, it is possible to take different Lifshitz scaling 
for space and time, and the resultant theory is called Horava-Lifshitz gravity [llQ- In the IR limit this gravity reduces 
to general relativity, and so it can be considered as a UV completion of general relativity. Motivated by this work on 
Horava-Lifshitz gravity, different Lifshitz scaling for space and time have been considered for type IIA string theory 
ll, type IIB string theory Q, AdS/CFT correspondence [1-Q, dilatonic black branes and dilatonic black holes 

111 [l^ . Another UV completion theory of general relativity which reduces to general relativity in the IR limit is 

called gravity’s rainbow [Q. In fact, it has been demonstrated that gravity’s rainbow is related to Horava-Lifshitz 
gravity 1^ . This is because both of these theories are based on the modifying the usual energy-momentum dispersion 
relation in the UV limit such that it reduces to the usual energy-momentum dispersion relation in the IR limit. It 
may be noted that such modification of the usual energy-momentum has also been obtained in discrete spacetime [l5l| , 
spacetime foa m ll^ . spin-network in loop quantum gravity (LQG) 53 , ghost condensation [l8| and non-commutative 
geometry [S |20||. The non-commutative geometry occurs due to background fluxes in string theory (^ . , and is 

used to derive one of the most important rainbow functions in gravity’s rainbow . 

It may be noted that the UV modification of the usual energy-momentum relation implies the breaking of the 
Lorentz symmetry in the UV limit of the theory. The spontaneous breaking of the Lorentz symmetry can occur in 
string theory because of the existence of an unstable perturbative string vacuum [^ . It is possible for a tachyon 
field to have wrong sign for its mass squared in string field theory, and this causes the perturbative string vacuum 
to become unstable. The theory becomes ill-defined if the vacuum expectation value of the tachyon field is infinite. 
It is also possible for the vacuum expectation value of the tachyon field to be finite and negative. In this case, the 
coefficient of the quadratic term for the massless vector field is nonzero and negative, and this breaks the Lorentz 
symmetry. The spontaneous breaking of the Lorentz symmetry in string theory has also been investigated using the 
gravitational version of the Higgs mechanism [2^. This has been done for the low energy effective action obtained 
from string theory. The Lorentz symmetry breaking has been also studied using black brane in type IIB string theory 
[?7| . In this analysis, moduli stabilization was studied using a KKLT-type moduli potential in the context of type 
IIB warped flux compactification. It was demonstrated that a Higgs phase for gravity will exist if all moduli are 
stabilized. Another study regarding the breaking of Lorentz symmetry was done by using compactification in string 
field theory [^ . It may be noted that various other approaches to quantum gravity also indicate that the Lorentz 
symmetry mig ht only be an effective symmetry which occurs in the IR limit of some fundamental theories of quantum 
gravity (29l433l |. 

Hence, there is a good motivation to study the UV deformation of geometries that occur in string theory. In fact. 
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motivate d by Lifshitz deformation of such geometries, and the relation between Horava-Lifshitz gravity and gravity’s 
rainbow |14j . recently rainbow deformation of geometries that occur in string theory has been performed. Thus, 
the modifications of the thermodynamics of black rings has been analyzed using gravity’s rainbow [s^ . It has been 
observed that a remnant exists for black rings in gravity’s rainbow. It has also been argued that a remnant might 
exist for all black objects in gravity’s rainbow [35|. This has been explicitly demonstrated for Kerr black holes, Kerr- 
Newman black holes in de Sitter space, charged AdS black holes, higher dimensional Kerr-AdS black holes and black 
Saturn (35j| . This was done by generalizing the work done on the thermodynamics of black holes in gravity’s rainbow 
[^ . As the usual uncertainty principle still holds in gravity’s rainbow [sj H^, and it is possible to obtain a lower 
bound on the energy E > l/Ax, using the usual uncertainty principle. This energy can be related to the energy of a 
particle emitted in Hawking radiation. Furthermore, the value of the uncertainty in position can be equated to the 
radius of the event horizon, E > I/Aa; « !/?’+■ This energy can be related to the energy at which spacetime is probed, 
and hence it describes the energy E in gravity’s rainbow. This is because effectively this particle emitted with the 
energy E can be viewed as a probe which is probing the geometry of the black hole. This consideration modifies the 
temperature of the black hole . The entropy and heat capacity of black hole in gravity’s rainbow can be calculated 
using this modified temperature. An interesting consequence of this modified solution is that it predicts the existence 
of remnants for the black hole. Thus, the temperature of the black hole reduces to zero, when black hole has a small 
but finite size. At this size, the black hole does not emit any Hawking radiation. The existence of black hole remnants 
can be used as a solution for the information paradox |^. Furthermore, it also solves problem related to the 
existence of naked singularity at the last stage of the evaporation of a black hole. As in this picture, a black hole does 
not evaporate completely producing a naked singularity, but rather a remnant is produced at the last stage of the 
evaporation of the black hole. The existence of a remnant also has phenomenological consequences. This is because 
it is not possible to produce black holes smaller than these remnants. This increases the energy at which mini black 
holes can be produced at the LHC [4l|. Recently, a lot of interest has been generated in gravity’s rainbow [42l - l^ . 
It may be noted that the rainbow functions have been constrained from experimental data Black hole solutions 
in gravity’s rainbow with nonlinear sources have been investigated in (d^ . In addition, the hydrostatic equilibrium 
equation for this gravity was obtained in Ref. [i^ . 

As there is a strong motivation to study rainbow deformation of geometries that occur in string theory, so we 
analyze the rainbow deformation of charged dilatonic black holes in this paper. It may be noted that dilaton gravity 
arises as a low-energy effective field theory of string theory [s^, IMl ■ The dilaton field is also a candidate for dark 
matter (s^ . In fact, in order to have better picture of nature of the dark energy, a new scalar field is added to the 
field content of the original theory . Black objects in presence of dilaton gravity, have also been investigated 

57l| . Recently, dilaton field has been used for analyzing compact objects and hydrostatic equilibrium of starts 
M. The evaporation of quantum black holes has also been investigated using two dimensional dilaton gravity 
6Ij . Motivated by these applications, we analyze dilaton field using the formalism of gravity’s rainbow. 

Thermodynamical aspects of black holes have been of a great interest ever since of pioneering works of Hawking 
and Beckenstein . The idea that geometrical aspects of black holes could be interpreted as thermodynamical 

quantities provides a deep insight into the connection between gravity and quantum mechanics. On the other hand, 
introduction and developments in gauge/gravity duality highlighted the importance of black holes thermodynamics 
In addition. Hawking and Page showed the existence of a phase transition for asymptotically anti de-Sitter 
black holes [t^. This phase transition was reconsidered through the use of AdS/CFT correspondence by Witten 
[80l |. These works motivated a large number of researches to be conducted in context of black holes thermodynamics, 
stability and their phase transitions [Sll - I^ . 

Recently, it has been demonstrated that, it is possible to treat the cosmological constant as the thermodynamic 
pressure in extended phase space. There are several reasons for such consideration which among them one can 
point out the existence of second order phase transitio n for black holes, Van der Waals like liquid/gas behavior in 
phase diagrams and formation of the triple point [IMoi. The consideration of the cosmological constant as a 
thermodynamical variable could be supported by studies that are conducted in context of AdS/CFT |l20l - ll23 |. In 
addition, it was shown that a case of ensemble dependency exists for charged 3- dime nsional black holes which could 
be removed by considering cosmological constant as a thermodynamical variable |l24j . The thermodynamical cr itical 



behavior of black holes in presence of different matter fields and gravities has been investigated in literature UMoi 


Another interesting method of studying thermodynamical structure of the black holes is through the use of geometry. 
It is proposed that one can build a phase space of the black holes by employing one of thermodynamical quantities 
of the black holes as thermodynamical potential and its corresponding extensive parameters as components of the 
phase space. The information regarding phase transitions of the black holes is within the singularities of Ricci scalar 
of the constructed phase space. In other words, the divergencies of the Ricci scalar of thermodynamical metric are 
representing bound a nd phase t ransition p oints . Th e thermodyna mical pot ential for this method could be mass which 


is used in Weinhold 


nd pnase t ransition p oints . In e tnermodyna mical pot 
fna [ml, Quevedo fTTl - lni| and HPEM fTil - lll7 | 


metrics or entropy which is employed in 


Ruppeiner jll8l , lll9l| metric. It was pointed out that Ruppeiner and Weinhold methods are related to each other with 
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temperature as conformal factor |ll2l - lll4 |. It was shown that for specific cases of black holes, Weinhold, Ruppeiner 
and Quevedo metrics may fail to pr ovide consistent results regarding phase transitions while the HPEM metric is 
proven to be successful one [11!tI - I117I |. In what follows, we use all of these methods to study phase transitions of the 
black holes. 

This paper is organized as follows. We obtain the charged black hole solutions in dilaton gravity’s rainbow and 
analyze their properties. This will be done by making the metric of charged black hole solutions in dilaton gravity 
depends on the energy. We also examine the first law of thermodynamics for this solution. Next, we study the stability 
of such solutions in gravity’s rainbow and phase transition of these black holes through heat capacity, geometrical 
thermodynamics and the analogy between cosmological constant and thermodynamical pressure. Finally, we obtain 
critical pressure and horizon radius through another method. Last section is devoted to conclusion. 


II. CHARGED DILATONIC BLACK HOLE SOLUTIONS IN GRAVITY’S RAINBOW 

In this section, we obtain charged black hole solutions in dilaton gravity’s rainbow and investigate their properties. 
This will be done by writing an energy dependent version of the metric for dilaton -Max well gravity. It may be 
noted that gravity’s rainbow is based on the generalization of doubly special relativity |l25l |. and so it is not possible 
for a particle to attain energy greater than the Planck energy in gravity’s rainbow. This is because gravity affects 
particles of different energies differently, and so the spacetime is represented by a family of energy dependent metrics 
in gravity’s rainbow [l^ . The gravity’s rainbow can be constructed by considering the following deformation of the 
standard energy-momentum relation 


E^f{e)-p^g\e)=m\ ( 1 ) 

where energy ratio is e = E/Ep^ in which E and Ep are, respectively, the energy of test particle and the Planck 
energy. The functions f{e) and g{e) are required to be constrained in such a way that the standard energy-momentum 
relation is obtained in the infrared limit. Thus, we require 

lini/(e) = l, lini 5 f(e) = l. (2) 

£->■0 E-S-0 

It may be noted that the spacetime is probed at the energy E, and by definition this cannot be greater than the 
Planck energy, /^(e) and g^{e) are called the rainbow functions and their functional forms ar e ph enomenologically 
motivated. Now it is possible to define an energy dependent deformation of the metric g{e) as |126l| 

g{e) = 77“'’ea(e) O eb{e), (3) 


where 


e«(,) = 

here cq and refer to the energy independent frame fields. 
The 4-dimensional action of charged dilaton gravity is 


5(e) 


127 | 


( 4 ) 


T = 


1 

IGtt 




n 


2(V$)^ 


U ($) - , 


( 5 ) 


where TZ is the Ricci scalar curvature, $ is the dilaton field and V ($) is a potential for $. The electromagnetic field 
is = di^A„ — in which A^ is the electromagnetic potential. In addition, it should be pointed out that a is a 
constant which determines the strength of coupling of the scalar and electromagnetic field. Due to the fact that we 
are looking for the black hole with a radial electric field {Ftr{r) = —Frt{r) ^ 0), the electromagnetic potential will be 
in the following form 


A^ = 5lh (r). 


( 6 ) 


Using variational principle and varying Eq. ([5]) with respect to the gravitational field g^^, the dilaton field $ and 
the gauge field we can obtain the following field equations 


Rau=2{ 


+2e-^“^ [F^r,F:i - , 


( 7 ) 
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( 8 ) 


= 0 . 


(9) 


In this paper we are attempting to obtain dilaton-Maxwell rainbow solutions. To do so, one can employ following 
static metric ansatz 




/ 2 (e) g‘^{e) 


dr^ 

'l'(r) 


+ r^R^ir)dnl 


( 10 ) 


where 'h(r) and R{r) are radial dependent functions which should be determined, and represents the line element 
of a 2— dimensional hypersurface with the constant curvature 2k and volume W 2 - We should note that the constant 
k indicates that the boundary of t = constant and r = constant can be a positive (elliptic), zero (flat) or negative 
(hyperbolic), constant curvature hypersurface with following explicit forms 


dnl = 


dO^ + sm Odip^, 
dO^ + sinh^ Od^p ^, 
dO'^ + dif^, 


k = 1 
k = -l 
fc = 0 


( 11 ) 


Using Eq. ®, one can obtain electromagnetic tensor as 


Ffr — 


qe 


2 a<J> 


^R{r) 


2 ’ 


( 12 ) 


where q is an integration constant which is related to the electric charge of the black hole. 

Here, in order to find consistent metric functions, we use a modified version of Liouville-type dilation potential with 
following form 


U($) 


2ka^ 


9^{e) 


+ 2Ae 


2 q$ 


(13) 


where ICij = i + ja^ and A is a free parameter which plays the role of the cosmological constant. It is worthwhile to 
mention that for the case of gie) = /(e) = 1 , one obtains 


lim U($) 

a{s)=f{e)^l 


2ka^ 2± 


-|- 2,A.€ 


2 a^» 


(14) 


which is the usua l Lio uville-type dilation potential that is used i n the context of Friedman-Robertson-Walker scalar 
field cosmologies [l28j| and Einstein-Maxwell-dilaton black holes [Ullili®- 

Next, we employ an ansatz, R{r) = in the field equations. The motivation for considering s uch a n ansatz is 

due to black string solutions of Einstein-Maxwell-dilaton gravity which was first introduced in Ref. [l3lj . Now, we 
are in a position to obtain metric functions. It is a matter of calculation to show that by using Eq. (113, the metric 
(fTUl) and the mentioned ansatz for R{r) , we have following solutions for the field equations, (Eqs. ([7]) and ([51)) 


4'(r) = 


/Ci.i 

/C_i4 



g^/Ci,i/^(£) /b 
52(e)/C_3_i \ry r2 


(15) 


= ( 10 ) 

where b is an arbitrary constant and 7 = a^/ICi.i- In the above expression, m is an integration constant which is 
related to the total mass of the black hole. It is notable that, in the absence of a non-trivial dilaton (a = 7 = 0), the 
solution (IT5|) reduces to 


m A 


J.2 ’ 


4'(r) = fc 


r 


3 g^{e) 


(17) 
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FIG. 1: ^'(r) versus r for fc = 1, m = 5, A = —0.5, b — 1.2 and q = 0.68. 

Left panel for a = 0.9, /^(e) = 1, = 0.85 (dashed line), g^{e) = 0.96 (bold line) and g^{e) = 1.20 (continuous line). 

Middle panel for a = 0.9, g^{e) = 1, f^{e) = 0.80 (dashed line), /^(e) = 1.05 (bold line) and /^(e) = 1.20 (continuous line). 
Right panel for g^{£) = 1.3, /^(e) = 1.3, a = 0.85 (dashed line), a = 0.87 (bold line) and a = 0.9, (continuous line). 


which describes a 4-dimensional asymptotically AdS topological charged black hole in gravity’s rainbow with a positive, 
zero or negative constant curvature hypersurface. 

In order to confirm black hole interpretation of the solutions, we look for the curvature singularity. To do so, we 
calculate the Kretschmann scalar. Calculations show that for finite values of radial coordinate, the Kretschmann 
scalar is finite. On the other hand, for very small and very large values of r, we obtain 


ViuY oc 

lim = 


_ 4 ^ 2,1 
^1,1 


12A(a4 - 2a2 + 2) fb 
^3,-1 U. 


47 


(18) 

(19) 


Equation dT51) confirms that there is an essential singularity located at r = 0, while Eq. (IT^ shows that for nonzero 
a, the asymptotical behavior of the solutions is not AdS. It is easy to show that the metric function may contain 
real positive roots (see Fig. [T|), and therefore, the curvature singularity can be covered with an event horizon and 
interpreted as a black hole. 


III. THERMODYNAMICAL QUANTITIES 


Now, we are in a position to calculate thermodynamic and conserved quantities of obtained solutions and examine 
the validity of the first law of thermodynamics. 

In order to obtain the temperature, we use the concept of surface gravity to show that the temperature of these 
solutions has following form 



(v) 

-27 

q^f{£) r+A 


1 ^ 1 

47 r 

rl 9 ^{£)f{£) 

U+J 

/{e)r+/C_i,i 


( 20 ) 


On the other hand, one can use the area law for extracting modified version of the entropy related to the Einsteinian 
class of black objects with the following structure 


'!^2rX / 
V(£) \r+) 


in which by setting a = 0 and g{£) = 1, the entropy of Einstein-Maxwell-dilation black holes is recovered. In order to 
find the total electric charge of the solutions, one can use the Gauss law. Calculating the flux of electric field helps 
us to find the total electric charge with the following form 

^ T^2qf{£) 

A-Kg{£) 


( 22 ) 
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Next, we are interested in obtaining the electric potential. Using following standard relation, one can obtain the 
electric potential at the event horizon with respect to the infinity as a reference 


U{r) = A^X^\ 


r —>-oo 




>r+ 


h. 

■ 


(23) 


Finally, according to the definition of mass due to Abbott and Deser 132l4l34l |. the total mass of the solution is 

W2h'^'^ 


M = 


87r/Ci,i5f(£)/(£) 


m. 


(24) 


It is worthwhile to mention th at fo r limiting case of g{e) = /(e) = 1 and a = 0, Eq. ([24|) reduces to mass of 
the Einstein-Maxwell black holes |l30l| . In addition, in obtained conserved and thermodynamical quantities, only the 
electric potential remains unaffected by considering gravity’s rainbow. 

Now, we are in a position to check the validity of the first law of thermodynamics. To do so, first, we calculate the 
geometrical mass, m, by using / (r = r_|_) = 0. Then by employing obtained relation for geometrical mass and Eq. 
for total mass of the black holes, we find 


M(r+,g) = 


W2A}A"< 


where 


A = 


-) 

r+/ 


27 


^1,-1 
^ 1,1 2 ^2 


Stt/Ci _ig(£)/(e)’ 


-^1,3 


^1,3 A ^3,1 / L \ 

— TCYJ I ^ \ 




g‘^{e)IC-3,i \r+J 
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It is a matter of calculation to show that 


=T & = 


\dS J 


\dQ J ^ 


= U. 


Therefore, we proved that the first law is valid as 

(dM\ 

" = ( 3S ) 


Q 




(25) 


(26) 


(27) 


IV. THERMAL STABILITY 


In this section, we study thermal stability of the solutions in context of the canonical ensemble. The stability 
conditions in the context of canonical ensemble are determined by the sign of the heat capacity. In other words, the 
positivity/negativity of the heat capacity is denoted as the black object being in stable/unstable state. Therefore, 
in order to study the stability of the charged black holes in dilatonic gravity’s rainbow, we study the changes in the 
sign of corresponding heat capacity. It is worthwhile to mention that investigating the behavior of the heat capacity 
enables one to obtain the phase transitions of the solutions at the same time. The root and divergence point of 
the heat capacity are denoted as bounded point (r+o) and second order phase transition point (r+c), respectively. 
Bounded point is related to the root of temperature and the sign of T changes at r+g, while we expect to obtain 
positive temperature at r+c. 

The system in canonical ensemble is considered to be in fixed charge. Therefore, we have 



-1 

Q 


(28) 


Considering the mentioned bounded and phase transition points, one can obtain 


r T = { 

dM\ _ Q 
dS JQ ~ ^ 

] (dT\ _ 

( d^M\ 

{ \dS)Q ~ 

\^)q ■ 


hounded point 
0 phase transition point 


( 29 ) 
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There are three valuable known cases for the rainbow functions which are characteristics of the rainbow solutions. 
These three cases are arisen from different phenomenological origins with an upper limit for considering the energy 
of test particle E 


E 

The first case is originated from loop quantum gravity and non-commutative geometry. In this case, we have 
following relations for rainbow functions of metric 

/(e) = 1, 5(e) = v^l-^e”. (30) 

The other case is constructed by considering the hard spectra from gamma-ray bursts which leads to [l^ 


- 1 

/(e) = —^, 5 (e) = 1. 

pe 


(31) 


Interestingly, opposite to previous case, in this one the effect of g{e) which is coupled to spatial coordinates of the 
metric is vanished, whereas in the case one which is related to loop quantum gravity the effect of coupling term for 
time component of the metric is vanished. 

Finally, in the case three, the choices of rainbow functions are due to constancy of the velocity of the light 


1 


/<"' = »(") = —A.' 

Using first law of thermodynamics, one can rewrite the relation for heat capacity into 


Cq 




dr 


dr+ J 


(32) 


(33) 


Q 


Now, by employing Eqs. (uni) and m with (j33p , on6 C3>ii show thcit hos^t CcipcLcity is 

(;t) ^ (^9^^-Fiff^(^)/^(£)+4^5 ^(£^)(ft) 


Cq = 


2g'^{e)IC-i^i (r^kg'^ie) - g'^IC3^ig'^{e)p{e) - 


47 


(34) 


Considering obtained relation for the heat capacity fEa. (1341) 1 and three mentioned cases for the rainbow functions 
of the metric (Eqs. (|301l32l) b we study the stability of the solutions. In case of horizon flat {k — 0), one can find that 
the root(s) of the heat capacity and divergence point (s) are given by following relations 


r+o = ^ - 


b^A 


^ 1,1 




(35) 


^ 1,1 

(36) 

Interestingly, for horizon flat, the root of the heat capacity, hence bound point is independent of the dilaton 
parameter, a while divergency of the heat capacity is a function of this parameter. On the other hand, in order 
to have positive real valued divergency in AdS spacetime, we have a > 1 restriction for dilaton parameter. This 
condition for dS spacetime is opposite. In other words, the real valued divergency is obtained if 0 < a < 1. (see Eqs. 
( 1551 ) and ( 1551 ) for more details). In order to charged black holes in dilatonic gravity’s rainbow with flat horizon be 
stable, following conditions must be hold 


f 64A/C-i,i \ 

V 5 ^(£)/^(£^) 5 ^^ 3 ,i/ 


I -Ar4/C_ni/Ci,i < q^gie)^ (;^) 








0.03 n 




FIG. 2: Cq (continuous line), T (dotted line) and TZ (dashed line) versus r+ for fc = 1, A = —1, fe = 5, i? = 1, Ep = 5 and 
q = 1. 

a = 0.9 (left panel) and a = 1.09 (right panel) ’’different scales”. 


As for the cases of A: = ±1, it was not possible to obtain the analytical relations for the root and divergence point 
of the heat capacity. Therefore, we employ numerical method for studying the properties of the heat capacity for 
spherical and hyperbolic horizons. As for the stability conditions, there are different orders of the horizon radius 
for each term. These terms will have dominant effect in specific regions of horizon radius and other parameters. 
Considering the effectiveness of these terms, the stability conditions will vary from one case to another one. Also, 
these effective behaviors may present different regions of stability and instability. Taking a closer look at different 
terms, one can see that the most effective parameter in stability conditions which modifies the exponent of horizon 
radius highly and changes the positivity and negativity of each term, is dilaton parameter, a. In other words, 
considering different values of a, stable/unstable regions will be modified highly. This highlights the effect of dilaton 
field on thermodynamical behavior of the solutions. 

In order to give a better insight regarding the thermodynamical behavior of these black holes, we study the behavior 
of the temperature. The reason is the fact that negativity of the temperature represents nonphysical systems which 
are not of our interest. Therefore, we study the conditions for positivity/negativity of the temperature. Considering 
Eq. (Uni), there are three terms which are related to electric charge, cosmological constant and topological structure 
of the metric. The effectiveness of each term is a function of their factors. Therefore, considering different values for 
these factors may lead to one of the following scenarios: one root and temperature is an increasing function of horizon 
radius (left panel of Fig. [2]), no root and temperature is negative with one maximum (right panel of Fig. , two 
roots with one region of positivity and two regions of negativity (Fig. [3|), one root and temperature is a decreasing 
function of the horizon radius (Fig. |4|). 

In general, charged term is always negative in Eq. (I20p . If one considers AdS solutions, the second term will be 
positive. As for the last term, if spherical solution is chosen, then in order for topological term be positive, the dilaton 
parameter must be a > 1. Whereas in case of hyperbolic horizon the condition will be modified into a < 1. On 
the other hand, if dS solutions are considered, the second term will be negative. Therefore, the possibility of having 
positive temperature depends on topological term with mentioned conditions for spherical and hyperbolic cases. It 
is worthwhile to mention that horizon flat of dS solutions has negative temperature. Therefore, it is not physical. 
Considering the mentioned changes in different terms of temperature, depending on the dominant regions of the 
different each term, the temperature will have positive/negative value with different behavior (which were pointed 
out in plotted diagrams). 

In order to elaborate the effects of the dilation field on thermal stability and mentioned behaviors for the temper¬ 
ature, we plot various diagrams using the first model of rainbow functions of gravity’s rainbow. 

It is evident that in case of temperature being an increasing function of horizon radius, for positive temperature, 
we have stable black holes and heat capacity is an increasing function of r_|_ (Fig. [5] left panel). Increasing the 
dilaton parameter will lead to formation of two stable and unstable states where both of these states are in negative 
temperature (Fig. [5] right panel). On the other hand, by increasing dilaton parameter, the temperature will have 
two regions of negativity and one positivity. In the positive region, a phase transition takes place between unstable 
larger state to smaller stable state. This phase transition point is represented by a divergency of the heat capacity. 
Increasing dilaton parameter leads to decreasing the place of divergency of the heat capacity and increasing the region 
in which we have unstable physical solutions (Fig. [3]). The larger root of the temperature is highly sensitive to 
variation of the a (see Fig. IS. For sufficiently large values of the dilaton parameter, the temperature is a decreasing 













FIG. 3: Cq (continuous line), T (dotted line) and TZ (dashed line) versus r+ for fc = 1, A = 
q = l. 

a = 1.1 (left panel) and a = 1.3 (middle and right panels) ’’different scales”. 




FIG. 4: Cq (continuous line), T (dotted line) and TZ (dashed line) versus r+ for fc = 1, A = 
q = l. 

a = 10 (left panel) and a = 15 (right panel) ’’different scales”. 
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FIG. 5: Cq (continuous line), T (dotted line) and TZ (dashed line) versus r+ for fc = 1, A = - 
and a = 1.1. 

Weinhold (left panel), Ruppeiner (middle panel) and Quevedo (middle panel) ’’different scales’ 
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function of the horizon radius with one root. The heat capacity in the region of the positive temperature is negative, 
therefore in this case rainbow solutions are unstable (Fig. 0]). It is worthwhile to mention that the root in this case 
is an increasing function of dilaton parameter (compare two diagrams of Fig. 0]). 


V. GEOMETRICAL THERMODYNAMICS 


In this section, we will conduct a study regarding the phase transition points of the obtained black holes through 
the use of geometrical thermodynamics concept. In this concept, one can construct the thermodynamical structure of 
the black hole through the use of thermodynamical variables. In other words, by using one of the thermodynamical 
variable as potential and its corresponding extensive parameters, it is possible to build phase space. The singularities 
of the Ricci scalar of this phase space, marks two different properties of the solutions: one is a bound point which is 
related to the root of the temperature and marks the physical and non-physical solutions. The other one is related to 
singularities of the heat capacity which are marking the points which system goes under phase transition. 

Considering such property for the constructed phase space, a valid approach of the geometrical thermodynamics 
produces a Ricci scalar which has singularities that cover both of the mentioned points. Depending on the thermody¬ 
namical potential, the extensive parameters would be different for each phase space. One of these potentials could be 
entropy which could be used to construct Ruppeiner phase space. An other potential could be mass which i s einploye d 
to build Weinhold, Quevedo and HPEM phase spaces. The mentioned phase spaces have following forms E g-lmli 



Weinhold 

-MTgYbdX^dX^ 

Ruppeiner 

{SMs + QMq) {-MssdS^ + MQQdQ^) 

Quevedo 

{-MssdS^ + MggdQ^) 

HPEM 


where their corresponding denominator of their Ricci scalars are 


(37) 


denom{TZ) = < 


{MssMqq - MIqY iS,Q) 
{MssMqq - MIqYt{S, Q)M^ ( 5 , Q) 
{SMs + QMQfMlsM^Q 


Weinhold 

Ruppeiner 

Quevedo 


HPEM 


(38) 


in which Mqq = ( 0 -)^, Msq = Mss = and Ms = {^) q- 

Now, by using Eqs. (E]), dm), (05]) and 02|), one can construct mentioned phase spaces and calculate their 
corresponding curvatnre scalar for these black holes. Due to economical reasons, we will not present obtained relation 
for Ricci scalar but demonstrate results in plotted diagrams (see Figs. 0101 for HPEM metric and Eig. 0] for other 
metrics). Eig. [5] shows that for specific values one can find cases in which Weinhold (left panel of Eig. [5|), Ruppeiner 
(middle panel of Eig. 0]) and Quevedo (right panel Eig. 0|) will not produce suitable divergencies in their Ricci scalar 
to cover mentioned points. In other words, the divergencies of their Ricci scalar may not coincide with root and 
divergencies of the heat capacity. On the other hand, it is seen that all the divergencies of the curvature scalar of 
the HPEM metric match with bound and phase transition points of the heat capacity. The nature of the behavior 
of the Ricci scalar around each o ne of thes e divergencies enables one to recognize whether it is a bound point or a 
divergence point in heat capacity [IIMil- 


VI. PHASE TRANSITIONS IN EXTENDED PHASE SPACE 

In this section, we investigate the existence of second order phase transition through the analogy between negative 
cosmological constant and thermodynamical pressure. The usual relation for pressure and cosmological constant is 
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give by (o^lOlL Il36l| 


P = - 


A 

Stt 


(39) 


It was shown that the gravitational theory under consideration may affect this relation and modifies it [l37L Il38j| . 
In calculations of the conserved and thermodynamical quantities, we found that these quantities were modified due 
to existence of gravity’s rainbow and dilaton field. It is natural to question whether the usual relation between 
cosmological constant and thermodynamical pressure could be modified in presence of the dilaton field as well as 
rainbow functions. To investigate such modification, we use the right hand side of the Eq. 0. It is a matter of 
calculation to show that (after removing parts related to electromagnetic field) 


t; oc A 



(40) 


Obtained relation indicates that although both dilaton field and rainbow functions modified thermodynamical 
quantities, only the dilatonic part has direct effect on the relation between cosmological constant and pressure. 
Therefore, we use following analogy for studying the critical behavior of the system 


A (JlP 

Stt \r+J 


(41) 


The conjugating quantity related to pressure is obtained through the use of enthalpy 


V = 



(42) 


Since consideration of the cosmological constant extends our thermodynamical phase space, the mass term plays 
the role of enthalpy. Therefore, by using Eqs. 65) and (1421) . one can find modified volume of these dilatonic 

black holes as 


V = 


/Ci,i 


^5,3 


(43) 


Clearly, the volume of these black holes is a function of both rainbow functions and dilaton parameter. In other 
words, contrary to some specific modified gravities, in this gravity, the volume of the black hole is affected by the 
presence of rainbow and dilaton gravities. Here, in order to have a positive and non-zero volume, we find a restriction 
—3 < a < 3. Since we are not interested in negative values of the a, we restrict ourselves to 0 < a < 3. 

It should be pointed out that due to the relation between volume of the black hole and horizon radius, one is able to 
introduce specific volume for these black holes which enables us to use horizon radius instead of volume in following 
calculations. So, the pressure is given by 


P = 


^3-i5(e)/(£) lC3-ig'^{e)[p{e)lC-i.iq^+r\] 


2 /Ci 1/C3 1 


Sr^Tr/C-ij/Ca,! 



(44) 


In order to find a relation for calculating critical volume, hence critical horizon radius, we use the concept of 
inflection point. In this method, one uses 




to find critical horizon radius which in case of this thermodynamical system is 

Tc = g/(e)\//C3,i/C2.i, 


which will lead to following critical temperature and pressure 

rj, ^ A:i,ig(£)re'"^-^ / b \ 

/3(e)g27r/Ci _i/Ci \rc) 


8^/2(£)g2/C3 • 


(45) 


(46) 

(47) 

(48) 
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FIG. 6 : P — r+ (left), T — r+ (middle) and G — T (right) diagrams for q = 0.1, h — 1, a = 0.7, g{e) = f{e) = 0.9 (continuous 
line), g{e) = f{e) = 1 (dotted line) and g{e) = /(e) = 1.1 (dashed line). 

P — r+ diagram for T = Tc, T — r+ diagram for P — Pc and G — T diagram for P — 0.5Pc- 


It is worthwhile to mention that the restriction that was observed was originated only from dilatonic part of the 
solutions. In other words, we have no restriction on values that charge and rainbow functions can acquire and our 
system is only thermodynamically restricted by dilaton parameter. 

Using obtained critical values, one can find the following critical ratio 


Pare 


gg(£)/^(£)^-i.i^3p 


\ 27 /Ci.-i 


(49) 


which shows that this critical ratio was modified due to the presence of dilaton held as well as rainbow functions. It 
is worthwhile to mention that critical horizon radius depends only on one of the rainbow functions whereas the other 
critical values and also the ratio are functions of both of them. In addition, it is notable that in the absence of 
dilaton held {a = 0) and low en ergy limit (/(e) = ^(e) = 1), Eq. (1491) reduces to the usual universal ratio in four 
dimensional Einstein gravity [HI, l99l | . 

Next, using the renewed role of the total mass of the black holes, we have Gibbs free energy as 


G = H-TS = M-TS. 

which by using Eqs. dini), dm, dm and (1411) will be 

-27 " 


G = 


ICl 1 


oi‘‘+3 

^1,1 


1C: 


2/C_3.i/(e)53(e) 


P 


w/C) (i)' 


^1,1 


6^7 


ii) 


-27 


^1,-1 

^1,1 


&27 


167rg(e) 


167r/(e)g(e) 


(50) 


(51) 


In order to see whether obtained critical values represent a second order phase transition, we study phase diagrams 
(P — r+, T — r+ and G — T diagrams) in Figs. IBUSl 

It is evident that for specihe values of different parameters, a second order phase transition is observed for obtained 
critical values (see Figs. |6]and[7l). The critical pressure (left panels of Figs. [6] and [7]), temperature and subcritical 
isobars (middle panels of Figs. in]and[71), energy of different phases and size of swallow-tails (right panels of Figs. [Bland 
0 are functions of gravity’s rainbow and dilaton parameter. The effects of rainbow functions and dilaton parameter 
on critical values are different from each other (compare Figs. [6] with 0. 

Interestingly, for a set of values, it is possible to obtain positive critical pressure and horizon radius whereas the 
temperature is negative. The plotted diagrams for these cases show a normal critical behavior in P — r+ diagram (left 
panel of Fig. 0 whereas in G — T diagram an abnormal behavior is observed (right panel of Fig. 0. In case of P — r+ 
diagram, also normal critical behavior is observed except that this behavior is located in negative temperature. 


VII. PHASE TRANSITION POINTS THROUGH HEAT CAPACITY 


In this section, we will obtain critical points through a method which was developed in Ref. Eol. In this method, the 
denominator of the heat capacity is employed to obtain an explicit relation for thermodynamical pressure. Obtained 
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FIG. 7: P — r+ (left), T — r+ (middle) and G — T (right) diagrams for 5 = 0.1, b = 1, g(e) = /(e) = 0.9, a = 0.75 (continuous 
line), a = 0.76 (dotted line) and a = 0.77 (dashed line). 

P — r+ diagram for T = Tc, T — r+ diagram for P = Pc and G — T diagram for P = 0.4Pc. 



FIG. 8: P — r+ (left), T — r+ (middle) and G — T (right) diagrams for q — I, b — 20, g{£) = /(e) = 2.5, a = 0.75 (continuous 
line), a = 0.76 (dotted line) and a = 0.77 (dashed line). 

P — r+ diagram ior T — Tc, T — r+ diagram for P — Pc and G — T diagram for P — 0.5Pc. 




FIG. 9: P versus r+ diagrams for q = 0.1, b — 1 and k — 1. 

left panel: g{e) = /(e) = 0.9 and a = 0.75 (bold continues line), P — 0.0455693 (continues line), a = 0.76 (bold dotted line), 
P = 0.0438404 (dotted line), a = 0.77 (bold dashed line), P — 0.0421757 (dashed line). 

middle panel: g{e) = /(e) = 0.9 and a — 1.65 (bold continues line), P — 0.0000080 (continues line), a = 1.66 (bold dotted 
line), P = 0.0000069 (dotted line), a = 1.67 (bold dashed line), P = 0.0000059 (dashed line). 

right panel: a — 0.7 and g{e) = /(e) = 0.9 (bold continues line), P = 0.0552557 (continues line), g{e) = /(e) = 1 (bold dotted 
line), P = 0.0592206 (dotted line), g{£) = /(e) = 1.1 (bold dashed line), P — 0.0630518 (dashed line). 
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relation may yield a maximum(s) for pressure which is(are) critical pressure(s) in which a second order phase transition 
takes place. This critical pressure is exactly the same as that of obtained through the use of phase diagrams. 

Using Eqs. (IM)) and (TO and solving the denominator with respect to thermodynamical pressure will lead to 
following explicit relation for pressure 

p _ ^-3,1 [gV^(g)/C3,i - rl] g^{e) 

Snr'llCs.ilCi^i 

It is evident that this relation is different from previously obtained relation for pressure (Eq. IHl) . Now, by using 
values that are employed for plotting phase diagrams (Figs. [ 5115 ]), we plot following diagrams (Fig. [ 51 ). A simple 
comparison shows that the maximums of the plotted diagrams are exactly where corresponding critical pressure and 
horizon radius are located in phase diagrams. This shows that these two approaches yield consisting picture regarding 
the critical behavior of these black holes. On the other hand, plotted diagram which corresponds to one with abnormal 
behavior (middle panel of Fig. also represents the characteristic behavior of the phase transition point. Therefore, 
in case of these black holes, a phase transition occurs in the mentioned critical point. 



VIII. CONCLUSION 

In this paper, we studied 4-dimensional charged dilatonic black holes in gravity’s rainbow and their thermal stability 
conditions. We obtained thermodynamical quantities such as temperature, electric charge, entropy and total mass of 
the black holes. These quantities were modified in gravity’s rainbow and became energy dependent. 

Next, we conducted a study regarding physical/nonphysical black holes (positivity/negativity of temperature) and 
thermal stability of the solutions. It was pointed out that dominant factor in studying these properties is the dilaton 
parameter. In other words, these properties were highly sensitive to variation of a. Due to different factors of dilaton 
parameter, different types of behavior were observed for temperature which put restrictions on the solutions being 
physical. Observed behaviors for the temperature were: a) two roots with maximum, b) an increasing (a decreasing) 
function of horizon radius with one root c) a negative definite function with one maximum. 

The analyzed behaviors were, increasing function of horizon radius with one root, increasing and decreasing function 
of horizon radius with two roots, decreasing function of the with one root located and being negative with one 
maximum located at negative temperature. 

As for the stability and phase transition, we found depending on the behavior of the temperature, heat capacity 
could have phase transition and stable state for larger values of horizon radius. In case of two roots for temperature, 
interestingly, a phase transition of larger/smaller black hole was observed. Finally, as for temperature being decreasing 
function of r+, for physical solutions (positive temperature), unstable solutions were observed. In other words, in this 
case, physical solutions are unstable. 

Next, geometrical approach was employed to study the bound and phase transition points of these black holes. 
It was demonstrated that the Ricci scalars of the phase spaces of Weinhold, Ruppeiner and Quevedo metrics, have 
divergencies which do not match with mentioned points while the singular points of curvature scalar of the HPEM 
coincide with roots and divergence points of the heat capacity. 

We also, studied the critical behavior of these black holes in extended phase space. It was shown that the usual 
relation between cosmological constant and thermodynamical pressure was modified due to existence of dilaton gravity 
whereas such modification was not seen for gravity’s rainbow. On the contrary, it was shown that volume of the black 
holes depends on both of these modifications. 

Then, we showed that the critical values are related to the rainbow functions as well as the dilaton parameter. In 
order to have positive critical temperature, we found restrictions which were purely dilatonic dependant. Therefore, 
one is not free to choose any value for dilaton parameter. 

In studying phase diagrams, two different behaviors were observed for different diagrams, especially in G — T 
diagrams. Interestingly, although we observed an anomaly in these diagrams, other corresponding phase diagrams 
presented usual thermodynamical behavior around critical points. In other words, the observed abnormal behaviors 
in phase diagrams present the existence of a second order pha se tra nsition for these black holes. 

Next, we used a new method which was introduced in Ref. [I04| . for studying the critical behavior of the system. 
This method is based on obtaining an explicit form for thermodynamical pressure from denominator of the heat 
capacity. It was seen that the maximum maximum of this relation is located at the critical pressure and horizon 
radius in which second order phase transition takes place. It was shown that in this method, for irregular behavior 
which was observed in phase diagrams, also a second order phase transition occurs. This indicates that these points 
are phase transition point despite their abnormal behavior. 

Finally, it is worthwhile to think about the physical interpretation of abnormal behavior which was seen in this pa¬ 
per. It is notable that, one can generalize obtained linear solutions in this paper to nonlinear case of electrodynamics 
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and investigate the effects of nonlinearity [m. In addition, one may investigate the extended phase space and thermo- 
dynamic criticality in higher order Lovelock-Maxwell gravity’s rainbow as well as Lovelock-nonlinear electrodynamics 
[l4nl Il4l| . These subjects are under examination. 
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